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Abstract 

This paper is devoted to the study of the behavior of the sohition ug e Hq{VI), as 
5 goes to 0, to the equation 

div(e5v4Vw) + k^ea^u = e^f in fi, 

where is a smooth open subset of with d = 2 or 3, / G L^(r2), fc is a non- negative 
constant, ^ is a uniformly eUiptic matrix function, E is a real function bounded above 
and below by positive constants, and eg is a complex function whose the real part 
takes the value 1 and —1, and the imaginary part is positive and converges to as (5 
goes to 0. Under some additional general assumptions on A and S, we characterize 
conditions on / under which ||u5||_f/i(n) remains bounded as 5 goes to 0. Under these 
conditions, we also show that us converges weakly in H^ijl) to a limit which is the 
solution to the limit equation; moreover, we obtain a formula for computing the limit. 
The applications of these results for perfect lens, cloaking, and illusion optics using 
negative index materials will be given. 

1 Introduction and the statements of the results 

1.1 Introduction and the geometric setting 

Negative index materials were first investigated theoretically by Veselago in |27| and 
were innovated by Pendry in |20) . The the existence of such materials was confirmed 
by Shelby, Smith, and Schultz in |25j (see also |26|). In [20j, Pendry showed that the 
negative index material slab (considered in Veselago's paper) acts as a lens not only for 
the propagating waves (for which the ray analysis of Veselago is valid) but also for the 
evanescent near-field radiation. In an earlier work [T7], Nicorovici, McPhedran, and Milton 
showed that negative index materials can be used for a (magnified) cylindrical perfect 
lens in 2d quasi-static regime. In [21], Pendry showed how a cylinder of material can be 
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designed to magnify an image with the same perfection through conformal transformations 
from the setting in [20j . In [24j , Anathana Ramakrishma and Pendry studied perfect lens 
solutions to spherical shells for spherical shells of negative index materials. In [llj, Lai et. 
al. investigate the way to create illusion optics by negative index materials again. In [lOj . 
Lai et. al. studied cloaking using negative index materials. More applications related to 
negative index materials can be found in references mentioned in the above works. 

In this paper, we extend several results in the works mentioned above to a more general 
setting for the wave equation in the time harmonic regime. More precisely, we study the 
behavior of the solution us G Hq{Q), as 5 goes to 0, to the equation 

diY{esAVu) + k'^soT.u = eof in 17, 

where 17 is a smooth open subset of M"' with d = 2 or 3, / G L^(r2), is a non- negative 
constant, ^ is a uniformly elliptic matrix function, E is a real function bounded above and 
below by positive constants, and es is a complex function whose the real part takes the 
value 1 and —1, and the imaginary part is positive and converges to as (5 goes to 0. Under 
some additional general assumptions on A and S, which are very related to the concept 
of complementary media introduced in (see also [lO]), we characterize conditions on 
/ G under which ||it(5||iyi(n) remains bounded as S goes to 0. Under these conditions, 

we also show that us converges weakly in H^{^) and strongly in L^(il) to a limit which 
is the solution to the limit equation; moreover, we obtain a formula for computing the 
limit. The applications of these results for perfect lens, cloaking, and illusion optics using 
negative index materials will be given. Our results are valid in the quasi-static regime 
(Theorem [T]) as well as in the finite frequency regime (Theorems [2] and Al ) for the 2d and 



3d cases. Our approach is based on transformation optics and the calculus of variations 
(or compactness arguments). Concerning transformation optics, we not only use the fact 
that the wave equation is invariant under changes of variables but also its Neumann 
boundary condition (see Proposition [l]) Our approach (see the proof of Theorems [l] 
and [2] in Sections 2.2 and 3.2) reflects the fact that negative index materials would make 



the geometry folded as mentioned e.g. in [13]. Our analysis could be generalized to the 
vectorial (Maxwell) case which would be considered elsewhere. Recently, in [2], the author 
characterize sources for which cloaking due to anomalous localized resonance appears in 
the 2d quasi-static case with A = I. The context in this paper and the one in [2] coincide 
only in the radial symmetric case. In that case, we emphasize that there is a difference in 
our characterization and the characterization given there. In this paper, we characterize 
sources for which the total energy of the fields remains bounded as the loss parameter goes 
to zero. In [2j, the authors characterize sources for which the dissipation energy blows up 
and the fields remains bounded outside a bounded region as the loss parameter goes to 
zero. We note that the dissipation energy goes to infinity implies the total energy goes 
to infinity; however the dissipation energy might be finite even though the total energy 
blows up (see [13] and Corollary [TJ. 



^In some specific situations, the normal derivative on the boundary changes its sign after a change of 



variables (see (1.25l). This fact will play an important role in our analysis. 
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Let us describe the problem more precisely. Let J7, f^i, and be smooth 

open subsets of M'^ {d = 2, 3) such that they are connected]^ 

0,1 CC CC 0,3 CC r24, CC ri, and Q is bounded. 

We assume that there exist xi £ l^i and X2 S f^s, and two diffeomorphisms 

F ■M2\{xi} ^n^\n2, G ■.n4\03^03\{x2} 

such that there exist diffeomorphism extensions of F and G on a neighborhood of dQ2 
and of di^3 respectively, 

F{Q2\Oi) = n3\ri2, 

F{x) = X on dil.2 and G{x) = x on dQ^, 

and 

G o F : r^i — 7- is a diffeomorphism if one sets G o F{xi) = X2- (1-1) 
The geometry of the problem is given in Figure 1 . 




Figure 1: Geometry of the setting 

In what follows, Bj. denotes the ball centered at with radius r. Here are some 
examples on the geometry. The first one is 

Example 1. Let m,l > 1, < ri < r2 < < r and let ^1 = Br, ili = Br^, il,2 = Br2, 
= Br^, 0,4 = M'^, and let yo £ ^^3- Define 

Fix) = , ^ X and Gix) = r^x. 



■^In this paper, for two subsets Di and D2 of R'', D\ CC D2 means Di C -D2, where D\ is the closure 
of Dx 
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Assume that = r^/r^ ^ . Then all conditions are verified for xi = X2 = 0. 

We next give the second example in which Q/j^ is bounded. 

Example 2. Let < ri < r2 < r^ < r^, r and let = Br, i^i = B^^, ^2 = -B^j, J^s = B^^, 
= Bj.^ ■ Define 

s / (^2 - ?"4)ri - (rs - r4)r2 I ,2 (r2 - r4)rf - (rs - r4)r| \x 

F{x)=[ 2 2 2 2 ^k|+r4h— T' 

V riTg — r^r2 rir2 — r^r2 / \x\ 

and 

G(x) = f-^^N+ ^^^^ 



r4 - rs r4 - rs/ |x| 
conditions are verified for xi = X2 = 0. 

For 5 > 0, £5 is defined as follows: 

f -l + i5 ifn2\f^i, 
esix) := ^ (1.2) 

[ 1 otherwise. 

Let A be a measurable matrix function and S be a measurable real function defined on Q, 
such that 

\\i?<{A{x%i)<m'' veeM", (1.3) 

for a.e. x G 17, for some < A < +00, and 

< ess inf S < ess sup S < +co. (1.4) 

^ n 

We are interested in the behavior of the solution us G H^{Q,) to the system 

dw{esAVus) + k'^£QT,us = eof in Q, 

Us = on dfl, 

as 6 goes to In this equation, the loss parameter involves only in the divergence term. 
The general case where the loss is also contained in the zero order term will be (briefly) 
discussed in Sections [LSI 



Our analysis in the zero frequency regime {k = 0) and in the finite frequency regime 
(k > 0) share the same approach although technicality is different. For the clarity, we 
state and prove the results separately. Before presenting these results in the next two 
subsections, let us introduce the following definitions: 

f A ifx£n\n3, . f s iix£n\n3, 

A:={ ^■.= { (1.5) 



■^We only consider the Dirichlet boundary condition. Other standard types of boundary conditions 
could be analyzed similarly. 
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and 

G^F^f iix e Q3. 
Here and in what follows, we use the standard notations: 



J^*A{y) = ^ ^ \ —, -F,S(y) = -f4, and TJ{y) 



J{x) 

where x = T^^{y) and J{x) = \ detDJ^(x)|. 



J(x)' 



J{x) 



(1.6) 



(1.7) 



1.2 Statements of the results in the quasi-static regime 

In this section, we are interested in the behavior of us, the unique solution in //q(J7) 
to the equation 

di\{e5AVu5) = eof in (1.8) 

as (5 is a small positive number. 



Remark 1. The existence and uniqueness of us will be established in Section 2.1. An a 
priori estimate of us will be also given there. 

Here is the main result of this section. 

Theorem 1. Let 6 > 0, f £ L^(r2) and us G Hq{^) be the unique solution to the equation 
(1.8). Assume A = F^:A in 0,^ \ Q2 o-nd in the case d = 3 assume in addition that 
A G Ci(Jl3 \f^2). We have 

a) Case 1: f is compatible with the system (see Definition ^below). Then the sequence 
{us) converges weakly in H^{Q,) and strongly in L^(r2) to uq the unique solution in 
Hq[^) to the equation 

div(eo^Vuo) = eof in Q, (1.9) 
as (5 —7- 0. Moreover uq = NI{f) where NI{f) is given by (1.13) 

b) Case 2: f is not compatible with the system. We have 



lim WusUhhu) = 
In the statement of Theorem [l| we use the following 



:i.io) 



Definition 1. Assume A = F^:A in Q-^ \ Q2 O'nd in the case d = 3 assume in addition 
that A G C^{Q,3 \ ^2)- We say that f G Lp'{Q) is compatible with the system if there exist 
U G H^in^ \ n2) and V G H^in^ \ ^2) such that 



dw{Avu) = F^f - f in n3\n2, 

U = on dn2, 

AVU-rj = ondn2, 



and < 



div(^Vy) = / 
V = W 



ext 



AW ■ ri = AVW ■ 7] 



ext 



in \ 0.2, 

on 51^3, 

on dO-i. 

(1.11) 
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Here W G H^{Q \ 80,-^) is the unique solution to the system 



div{AVW) = f 
W = 
[W] = -u 

[AVW • rj] = -AVU ■ r? 



in Q\ d^ls , 
on di} , 
on 80,3 , 
on 80,3 , 



(1.12) 



If f is compatible with the system then NI{f) is defined as follows 

w ifx€n\n3, 



NI{f) 



V 

{U + V)oF 
W oGoF 



if X &VL3\ 9,2, 

if X £ \ ^1, 
if X £ Qi. 



(1.13) 



We recall here that A and / are given in ( 1.5 ) and ( 1.6 ). In this paper, [ip] := ipcxt — fmt 
on 8D, and rj denotes an outward normal unit vector on 8D for any smooth open bounded 
subset L> of M°'. 

Remark 2. The uniqueness of U and V follow from the unique continuation principle 
(see e.g. \22^ and JJ^) and the fact that A G 0^(0,3 \ 0,2) in the 3d case. The existence 
and uniqueness of W are standard and can be proved by Fredholm's theory. The details 
are left to the reader. 

Remark 3. In general, ||'iio||_H'i(o) is not uniformly bounded with respect to ||/||l2 i.e., 
there does not exist a positive constant C , depending only on the geometry, S, an A, such 
that 

\\uo\\m<C\\f\\L2. 

This is a consequence of the fact that \\us\\ui is bounded only for the class of compatible 
functions f G L'^{Q) which is not open in L'^{Cl) in general. Hence there is a difference on 
the stability between equation (1.8) with 5 = and standard elliptic equations of divergence 
form. 

Remark 4. We note that since = ^ in \ 0,2, U + V {£ H^{Q3 \ ^2) satisfies: 

F^f in Q3\0,2, 
on 80,2, 
on 80,2. 



div [F^AV{U + V)^ 

u + v = v 

{ F^AV{U + V)-r] = F^AVV ■ ri 



From Theorem [T| it is clear that if / is compatible and / = in ^3, then uq = uqoGoF 
on 80,1 since C/ = 0. This is a phenomena of perfect lenses. More application of Theorem[T] 
on perfect lenses, cloaking, and illusion optics will be given in Section |4j The fact that 
in the context of Theorem [T| there are sources such that the energy of the fields goes to 
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infinity has been observed in several situations see e.g. [13], [6]. The media —A in 1^2 and 
F^:A in \ $^2 might be called complementary medium as suggested in ^21j (see also [10] 
and Section |4]). 

Assume f = in ft-^. From the definition of U ( |1.11[ ) and the unique continuation 
principle, it follows that U = in \ Q2. Therefore, by ( |1.12[ ), W e H^{i^) is the 
solution to the equation 

div(iVVF) = / in n. 

Assume in addition that A = ^ in \ O2 i.e. G^:F^:A = A in ^3 \ r22. It follows from 



the definition of V (1.11) and the unique continuation principle that V = W in fl3\fl2. 



Hence we obtain the following corollary 

Corollary 1. Let d = 2, 3 and f G L^(r2) be such that f = in Q^. Assume that G*-F*A = 
F^:A = A in Q,3\Cl2, and in the case d = 3 assume in addition that A £ C^{Cl3\Q2)- Then 
f is compatible and (us) C Hq{Q,), the solutions to the equation (1.8), converges weakly 



in H^{fl) and strongly in L'^{Q) to uq the unique solution in Hq{Q,) to the equation 

div(eo^V'Uo) = £0/ in ^) 

as (5 — )• 0. Moreover, 

Uq = W in 

where W G Hq{Q) is the unique solution to the equation 

div(iVVF) = f inn. 



We recall here that A is defined in (1.5). 



1.3 Statements of the results in the finite frequency regime 

In this section, we consider the finite frequency regime. We are interested in the 
behavior of us, the unique solution in Hq{Q) to the equation 

div{£sAVus) + k'^eoT.us = £of in Q, (1-14) 
as 5 is a small positive number. We assume in addition that the following systems 

div(^V';;) + k'^T.v = in \ 1^2 , 

v = ondn , (1.15) 

AVv ■ rj = on 80,2 , 

and 

f div(iV[/) + k'^tU = inn , 

(1.16) 

[ U = ondn, 

have only zero solution in H^{n \ ^12) and H^(n) respectively]^ 

''if we consider scattered waves, these conditions are not necessary. However, we need to impose the 
outgoing condition, see Appendix A. 
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Remark 5. The existence and uniqueness of us will be established in Section 3.1, An a 
priori estimate of us will be also given there. 

The main result in this section is 

Theorem 2. Let d = 2, 3, 5 > 0, / G L^(r2) and in the 3d case assume in addition that 
A G C^{0,3 \ 172)- Let Us G Hq{CI) be the unique solution to the equation (1.14). Assume 
that F^:A = A and F^:Y, = T, in \ 172- We have 

a) Case 1: f is compatible with the system. Then the sequence (ug) converges weakly 
in H^{Q) and strongly in L^(r2) to uq G Hq{Q), the unique solution to the equation 

div(eoVtio) + k'^eoT.UQ = eof in Cl. (1-17) 

as 5 — )• 0. Moreover, uq = NL(k, f) where NL{k, f) is given in (1.22) 

b) Case 2: f is not compatible with the system. We have 



lim \\us\\m(Q,) = +00. 



8^0 



:i.i8) 



In Theorem [2] we use the following 

Definition 2. We say that f G L^(il) is compatible with the system if there exist U G 
H^{n3 \ and V G H^i^s \ ^2) such that 

' div{AVU) + k^i:U = - f inQ3\n2, 

U = on dVt2, (1.19) 

AVU-ri = Q ondQ2, 



and 



( dw{AVV) + k'^YlV = f 

V = w 



cxt 



AW ■ T] = AVW ■ T] 



cxt 



in Vt-i \ CI2, 
on (903, 
on 



:i.2o) 



Here W G H^{Q \ dil^) is the unique solution to the system 



div(iVW^) + k'^tW = f 
W = Q 



[W] 
[AVW ■ ry] 



Then NI{k, f) is defined as follows 



NlikJ) 



-U 

-AVU ■ r] 

W 
V 



in 0,\ d^ls , 
on do, , 
on dfl^ , 
on . 

if X £ ft \ Qs, 
i/x G 1^3 \ 1^2, 



:i.21) 



{u + v)oF if xen2\ni, 

WoGoF ifxe^i. 



(1.22) 
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We recall here that A and / are given in (1.5) and (1.6). 



Remark 6. The uniqueness of U and V follow from the unique continuation principle 
see e.g. [2^/ and and the fact that A G C^(J]3 \ Q2) in the 3d case. The existence and 
uniqueness of W are standard and can be proved by Fredholm's theory since the system 
(1.16) is well-posed. 



Remark 7. We will present the results on scattered waves in the appendix. 



L2 



(see 



Remark 8. In general, ||wo||_ffi(n) '^ot uniformly bounded with respect to 
Remark^. 

Remark 9. We note that since F^A = A and F*S = in Q.z\^2, U + V e H^{^z \ ^2) 
and satisfies: 

' div (f^AV{\J + y)) + k^FJ^{U + = F,/ in ^3 \ ^2, 

U + V = V on 8^2, 

F^AV{U + V)-r] = F^AVV • rj on 8^2- 

Applications of Theorem [2] for perfect lenses, cloaking, and illusion optics will be 
given in Section |4j The media —A, — S in VI2 and F^A,F^T; in \ Q.2 might be called 
complementary medium as suggested in |21j (see also [10| and Section [4]). 

Assume / = in $^3. From the definition of C/ (1.19) and the unique continuation 
principle, it follows that f/ = in \ il2- Therefore, by ( |1.12[ ), W G -f^o(^^) is the 
solution to the equation 

div(iVVF) + k^tW = f in ft. 

Assume in addition that A = A and S = S in 1^3 \ 0,2 i-e. G*-F*A = A and = S in 

^3 \fl2. It follows from the definition of y (1.11) and the unique continuation principle 
that V = W mQ3\0,2- Hence we obtain the following corollary 

Corollary 2. Assume G^,F^A = F*^ = A in \ $^2; f = in ^^3, and in the 3d case 
assume in addition that A G C^{D,3 \ O2). Then f is compatible and (us), the solution to 
the equation (1.14) in Hq{Q), converges weakly in H^{Q,) and strongly in L'^{Q,) to uq the 
unique solution in Hq{^1) to the equation 

div(eo^Vtio) + /c^eoSuo = eo/ in 

as 5 — )• 0. Moreover 

Uq = W in 

where W G F[q{Q) is the unique solution to the equation 

div{AvW) + k'^tW = f inn. 



We recall here that A and S are defined in (1.5). 

In Theorem [2j we just consider the case in which the loss parameter involves only in 
the divergence term. In fact, the analysis in this paper can be extended to the general 
case in which the loss parameter involves in both terms. We have 
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Theorem 3. Let d = 2, 3, 5 > 0, / G L^(r2) and in the 3d case assume in addition that 
A G C^{Q3 \ f^2)- Let Us G Hq{Q) be the unique solution to the equation 

dw{esA\7us) + k^esT.us = f. 

Assume that F^^A = A and -F*S = S in f^s \ We have 

a) Case 1: f is compatible with the system. Then the sequence (us) converges weakly 
in H^{U) and strongly in L^(r2) to uq G Hq{Q), the unique solution to the equation 

div(eoVuo) + k'^eoT.UQ = eof in 0. (1.23) 



as 5 — )• 0. Moreover, uq = NL{k, f) where NL{k, f) is given in (1.22) 
b) Case 2: f is not compatible with the system. We have 

lim \\us\\H^n) = (1.24) 

The proof of Theorem |3] follows similarly as the one of Theorem [2] (see Section [s]). The 
details are left to the reader. 

1.4 Ideas of the proofs of Theorems [l] and [2] 

In this section we describe briefly the ideas of the proofs of Theorems [l] and [2] The 
proof of Theorem [T] is proceeded as follows. We first show that if / is compatible then 
u G Hq{Q,) is a solution to the equation ( |1.8[ ) if and only if u = NI{f). This is the 
crucial step of the proof. The key ingredients for the proof of this fact are properties of 
transformation optics stated in Proposition [T] below where we essentially use (1.25). Using 



the fact that NL{f) is the unique solution in Hq{Q) to the equation (1.8) with 5 = 0, we 
can easily establish the boundedness of us in H^{Q) by applying an elementary a priori 
estimate for the solution of equations of class (1.8) (Lemma [T| to the difference us — uq. 



We emphasize here that the constant in the estimate of Lemma [T] blows up as 1/5 as 5 goes 
to however, the source in the corresponding equation is of the order 5: it is exactly the 
place where we make use the existence of uq. The convergence of us follows immediately 
from its boundedness and the existence and uniqueness of uq . The proof of Theorem [2] 
follows in the same way. The only difference between these two proofs lies mainly on the 



way we establish a priori estimates for the solutions to the corresponding classes ((1.8) 



vs (1.14)); the proof in the finite frequency case is more complicated (see the proof of 
Propositions [2] and [s]) . This situation happens naturally and quite often since obtaining 
estimates for the finite frequency case usually requires more work than the quasi-static 
case. 

Remark 10. In this paper, we only consider structures generated from the class of trans- 
formations F and G mentioned in Section \l.l\ As a matter of fact, negative index material 
regions discussed in this paper are finite. Nevertheless, our analysis is quite robust and 
could be generalized for other situations e.g. the ones considered in JWj fl^ . and fB^. 
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We recall the following well-known proposition, which plays a very important role in 
our approach, on the transformation optics or the change of variables formula for the time 
harmonic wave equation. 

Proposition 1. Let k > 0, Di and D2 are two open subsets ofM.'^, T be a diffeomorphism 
fromDi ontoD2, a G [L°°(Di)]'^^'^ anda € L°°(i:>i). Fixu G H^{Di) andsetv = uoT-^. 
We have 

div(aVn) + k'^au = f 



iff 



dw{naVv) + k^nav = T*/ 



Assume that Ti and T2 are open subsets of dDi and dD2 such that Ti and T2 are smooth, 



T2 = T{Vi), and T := T Ti 



iff 

where 



r2 is a diffeomorphism. We have 
aVu ■ rji = gi on Ti 

%aVv ■ rj2 = g2 on T2, 
g2{y)=gi{x)/\detVT{x)\ withx = T~\y). 



Here r/i and r]2 are normal unit vector on Ti and T2 directed to the exterior of Di and 
D2. In particular, ifTi = T2, T(x) = x on Ti, D2n Di = 0. We have 



T^aVv ■ rji = —aVu ■ rji on Ti = T2- 



:i.25) 



Proof. Proposition [T] is a consequence of the change of variables. The first equivalent 
relation is well-known and can be proved by using the weak formula. The second equivalent 
follows similarly. A close variance and the approach mentioned here can be found in 
Section 2.2]. The details are left to the reader. □ 



Remark 11. Identity (1.25) gives a connection between negative index material and trans- 
formation optics. This identity will play an important role in our analysis. 



1.5 The organization of the paper 

The paper is organized as follows. In Section [2] and |3j we respectively establish Theo- 
rems [T] and [2j Section [4] is devoted to applications of Theorem [T] and [2] The results related 
to the scattered waves will be presented in the appendix. 

^In the identity below, VT stands for the gradient of a transformation from a (d — l)-manifold into a 
(d — l)-manifold, and det VT denotes the determinant of (d — 1) x (d — 1) matrix. 
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2 Quasi-static regime - Proof of Theorem [l] 

This section which contains two subsections is devoted to the proof of Theorem [T| In 
the first subsection, we estabhsh basis properties of the solution to equation ( |1.8| ) such 
as existence, uniqueness, and an a priori estimate. In the second subsection, we give the 
proof of Theorem of [T] 

2.1 Preliminaries 

The main result of this section is the following proposition which implies the existence 
and uniqueness of us, and give an a priori estimate for it. 

Proposition 2. Let d = 2,3,ge H-^{n) (the duality ofH^{n)), and A G [L~(J7)]^^'^. 
Assume that A satisfies (1.3). There exists a unique solution vs G Hq{Q) to the equation 

div{esAVvs) = g in ft. 

Moreover, we have 



IHWrnin) < Cmax{l,l/6}\\g\\H-i(^ay 
In this section, C denotes a positive constant depending only on and A. 
Proof. By Lax-Milgram's theorem, it suffices to verify that 



(2.1) 



/ es{AVu,Vu)dx >Cmin{5,l} / iVn^dx Vn G ^0^(0) (2.2) 
Jn Ju 



for some positive constant C. Indeed, we have 



es{AVu,Vu)dx >{a-bY + 6H^, 



where 



o = / {AVu, Vu) dx and & = / (AVu, Vn) dx. 

On the other hand, for 5 < 1/2 

(a - bf + S^b"^ =(1 - 5V2)a^ - 2ab + b'^/{l - 5^/2) + 6^a^/2 + {6^2 - <5V2)6V(1 - ^^2) 



> 



(1 - 6^/2)'/^a - 6/(1 - 6^/2)'/^] + C6\a^ + b^), 



and for 5 > 1/2, 



(a - bf + SH"^ > (a - bf + > C{a^ + 6^). 
Therefore, we obtain (|2.2|). 



□ 
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2.2 Proof of Theorem [T] 

Case 1: / is compatible. The proof is derived from the following steps: 
Step 1: Let v G Hq{Q) be a solution to the equation 

div(eo^Vu) = sof. 



(2.3) 



We prove that v = NL{f). 

Step 2: Define uq := NI{f). We prove that uq G Hq{Q) is a solution to the equation 

div(eo^Vuo) = eo/- 
Step 3: We prove that {us)o<:s<i is bounded in H^{^1). 

Step 4: We prove that (us) converges weakly in H^{il) and strongly in L^(il) to uq as 
6 goes to 0. 

We note immediately here that as a consequence of Steps 1 and 2, NI(f) is the unique 



solution in Hq{Q,) to the equation (1.8). We now process the proof of these steps. 
Step 1: Assume v G Hq{Q) is a solution to the equation 

div(eoAVw) = eo/ in ^■ 

Set 



vi = V o F ^ in $^4 \ 



and 



1 if G 1^3 \ i72, 
—1 if x G \ $73. 
From Proposition [ij vi G H^{Q3 \Q,2) <^ Hj-^^{fl4 \ 0,2) and vi satisfies 

div{iQF^A'Vvi) = e^F^f in $^4 \ ^2, 

and 

vi = V on dVt,2 and F^^AVvi ■ rj = AVv.rj 



(2.4) 
(2.5) 

(2.6) 



on 



cxt 



In the last identity, we use the fact that F^:A\7vi-r] = —A\7v-r] on 80,2 by Proposition 1 

int U 



and AVv ■ rj 



ext 



AVv ■ T] on 8^2 by the transmission condition on dVt2- Define 



U = fi — u in O3 \ ^2. 
Since F^A = A m. ^■i\Vt2-, it follows that 

' div(AVU) = F^f-f in J73 \ ^2, 



(2.7) 



U = 
AW • ?? = 



on dn2, 
on dQ2- 



(2. 
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Since A G C^{^3 \ ^2) in the 3d case, applying the unique continuation principle, see e.g., 
\22\ and , we have 

V = U inQ^X ^2. (2.9) 



We recall here that U is defined in (1.11). Define the function V2 on as follows 

[ -ui o G^^ if X G Q3, 
[ V iixe9.\^3. 

Using (1.1) and applying Proposition [T| we have 

dw{AVv2) = f m n\dQ3, 

and, on SOs, 



AVv2 • r] 



AVv2 • f] 



int 



It follows from (|2^ that 

AVv2 ■ r] 



ext 



AVv ■ 1] 
AVv ■ 1] 
AVv ■ 7] 

A'Vv2 ■ f] 



+ F^AVvi ■ r] 



oxt 



oxt 



F^AVvi ■ T] 



ext 



int 



(by (2.6)) 



F^AV{v + U) • r/ 



-F^AVU ■ ri 



int 



(by ^ 



Since G{x) = x on dil, we also obtain, on d^ls, 



V2 


- V2 


= V 


- Vl 


= V 


- Vl 


= V 




ext 


int 


oxt 


ext 


oxt 


int 



ext 



-U 



int 



A combination of ( |2.11[ ), ( |2.12| ), ( |2.13| ), and ( |1.12[ ) yields 

V2 = W in Q. 

This implies that v satisfies 

( dw{AVv) = f 
v = W 

A\/v ■ T] 



ext 



int 



AVW-T] 



oxt 



in Q3 \ O2, 
on 80,3, 
on 9^3. 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



Since A G 0^(03 \ 02) in the 3d case, applying the unique continuation principle, see e.g., 
[22] and [1], we have 

v = VmQ3\ri2. (2.15) 



14 



From (|2.10D, K7\, and it follows that 



V = < 



V2 in n \ O3, 

U + i;i inr23\r22, 

vio F in 02 \ f^i) 

V2 o G o F in 



We deduce from (2.9), (2.14), and (2.15) that 



V = < 



w mn\n3, 

V in \ 1^2, 

{V + U)oF in02\Sli, 

WoGoF in ill. 



Therefore, 



V = NL{f) in n. 

The proof of Step 1 completes. 
Step 2: We first prove that 

div(ylVuo) = / in (f] \ (li) \ {8^3 U dn2) 

where uq := NI{f). Since 



div(AvvF) = / in \ 

dW{AVV) = / in f]3 \ f]2, 



and 
it follows that 

div(^Vno) = fin (n\ 0.3^ U (^3 \ ^2 
Since = ^4 in 1^3 \ 1^2, U + V e H^{^l3 \ 0,2) and satisfies: 

div (^F^AV{U + V)^ = F^f in ^3 \ ^2 

we obtain 

div(^Vno) = / in 1^2 \ Oi. 



A combination of (2.17) and (2.18) yields (2.16). 



and 



It remains to verify 

[AVuo] = [uq] = on 8^3, [eoAVuo] = [uq] = on 80,2 

div{eoAVuo) = eof in $72- 



(2.16) 



(2.17) 



(2.18) 

(2.19) 
(2.20) 
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From the definition of V pTTTj ), of W and of NI{f) ([TlS]), we have 



[AVuq] = [no] = on 
Since C/ = and A^U • = on dVt2^ it follows from Proposition [T] that 

[eo^Vuo] = [uq] = on 80.2, 

Define 

Vi = 



V + U inl73\f^2. 



By Proposition [T| we have 

diviioF^VVi) = ioFJ in O4 \ ^^2, 
where eo is given in ( |2.5| ). From Proposition [T| it follows that 

div(eo^Vuo) = eo/ in O2. 

Step 3: Set 

vs = Us — uq in Q. 

We have, in fi, 

di-v{esAVvs) = diY{esAVus) — div(e5ylVMo) = £0/ — div(eo^Vno) + div [{eq — es)AVuo] . 
Since div(eoVuo) = Sof in we obtain 

div{esVvs) = div [(eo — £s)AVuq] in Q. 
We deduce from Proposition [2] that 

l|Vt;5||L2(o) < C||Vno||L2(n), 
which yields, since us = vs + uq, 

WVush^in) < CWVuoh^n)- 
Since us G Hq{Q), by Poincare's inequality, it follows that 

hsWmiU) < C'l|Vno||L2(s^). 

Step 3 completes. 

Step 4: The conclusion of Step 4 follows from Step 3 and the fact that the limit of us 



up to a subsequence) satisfies (2.3) and this equation has a unique solution in Hq{Q). 



Case 2: / is not compatible with the system. We process the proof by contradiction as 



follows. Assume that (1.10). Without loss of generality, there exists a sequence (us) (up 
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to an extraction) bounded in H^{Q) such that us converges weakly to u in H^{Q). It 
fohows that u G Hoi^) is a solution to the equation 

div(eo^Vti) = eo/ in ^■ 

Define 

U = uq o — uq in \ $^2 and V = uq in Qs\^2- 

This implies, as in Step 1, U and V satisfy ( |1.11[ ) respectively. We have a contradiction 
since / is not compatible with the system. □ 

Remark 12. Our analysis (mainly in Steps 1 and 2) would make media whose some 
region contains negative index materials behave like a folded space as noted in \14^ - 

3 Finite frequency regime - Proof of Theorem [2] 

This section which contains two subsections is devoted to the proof of Theorem [2] In 



the first subsection, we establish basis properties of the solution to equation ( 1.14[ ) such 
as existence, uniqueness, and an a priori estimate. In the second subsection, we give the 
proof of Theorem of[2j In this section, we always assume that (1.15) and (1.16) have only 
zero solution. 

3.1 Preliminaries 



In this section, we prove the the existence and uniqueness of the system (1.14) and 
establish an a priori estimate for its solutions. For this end, we need the following three 
lemmas. The first one is 

Lemma 1. Let d> 2, O be a smooth bounded open subset ofW^ and T C O be such that 
n'^-^ir) > (fl Assume u G H^{0). We have 

||«||l2(o) < c(^\\Vu\\l2^o) + lkllL2(r)), 
for some positive constant C independent of u. 

Proof. The proof follows from standard compactness arguments. The details are left to 
the reader. □ 

The second one is from [TJ Lemma A.l]: 

Lemma 2. Let d > 2, O be a smooth bounded open subset o/M'^, and u G Hdiv{0), i.e., 
u G [L'^{0)Y and divu G L'^{0). We have 

\W ■ "11^-1/2(50) - C'll^lli2(o) (ll divu||i2(o) + ||^i||L2(o)) > 

for some positive constant C depending only on O. 
Sy^d-i (jgjjQ^gg ^}jg g( _ 1-Hausdorf measure. 
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Remark 13. The standard estimate in the i?div(0) is: 

\\u ■ n\\jj^i/2(^go) < C(^\\ divM||i2(o) + \\u\\l2{o))^ 

which is weaker than the estimate stated in Lemma and not good enough for the proof 
of Proposition^ Lemma^is also useful in other situations see e.g. ^ and jTSf . 

The third lemma is a consequence of Fredholm's theorem (see e.g., [H [T^]) and the 



fact that the system (1.15) has only zero solution. 



Lemma 3. Assume that (1.15) has only zero solution. Let g £ H ^^'^(8^12) o-nd v E 



H^(Q, \ 1^2) be the unique solution to the system 

div(^V?;) + k'^av = in (7 \ 1^2 



v = 
AVv -rj = 



on dQ , 
on ■ 



We have 



\\v\\min\n2) ^ C\\g\\H-i/2(^9Q^), 
for some positive constant C independent of g. 

We are ready to state the main result of this section which implies the existence and 
uniqueness of ug in Theorem [2] and provides an a priori estimate for it. 

Proposition 3. Let g G H^^{^}) (the duality of Hq{Q,)). There exists a unique solution 
vs G Hq{Q) to the equation 

dW{esA'Vvs) + k'^Eo'Evs = g in Q. 

Moreover, we have 

hsWmin) < c(^-\\g\\H-^n) + Wgh^ni) + ll9llL2(n2\f^i))' (^-l) 
for some positive constant C independent of g and 5, as 6 is small^ 

Remark 14. We note here that the conditions g G H~'^{^), g G L'^[^2 \ ^i); cind g G 
-L^(ili) do not imply g G L'^(Q2)- Here is an example. Let v G //g(r2) such that v G 
H'^{ni),v £ H'^{9.2\^i) hutv ^ H'^{VL2). Setg = Av. Then g e {Vt) , g e L'^ {^2\^i) , 
and g G L^(rii), hut g L'^{Q.2). Hence one cannot write \\g\\L'^{yi^^) + 115111,2(^22X^1) ^■^ 



See Remark 14 for comments on ( |3.1| ) 
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Proof. The existence of vs will follow from its uniqueness by Fredholm's theorem (see 
e.g. [H m]). We now establish the uniqueness by showing that vs = if vs € Hq{Q) is a 
solution to the equation 

div{esA\7vs) + k'^eo'Svs = in Q. 
Multiplying this equation by vs and integrating the obtained expression on Q, we have 

£s{AVvs,Vvs) - k'^eo^vsl'^ = 0. 

This implies, by considering the imaginary part, 

{AVvs,Vvs) = 0. 

It follows that vs is constant in \ ^i- Thus f 5 = in 1^2 \ ^^i since = div{£sA'Vvs) + 
k'^Eo'Evs = k'^Eo'Evs in \ f^i- As a consequence, 

div{AVvs) + k'^Evs = in il. 



This implies vs = in Q since (1.15) has only zero solution in H^{Q). The proof of the 
existence and uniqueness of vs completes. 



We now establish (3.1) by contradiction arguments. Assume that there exists (gs) C 
H-^{n) such that 

\H\\m{n) = 1 and ^H^^Hh-i + Wdsh^ini) + Wash^QiXni) ^ 0> (3-2) 
as 6 goes to 0. Here vs € Hq{Q,) is the unique solution to the equation 

div{esAVvs) + k'^eoEvs = gs in fl. (3.3) 
Multiplying this equation by vs and integrating the obtained expression on $7, we have 

/ es{AVvs,Vvs) - k'^eoE\vs\'^ = - '■ 
Jn Jn 

Considering the imaginary part and using the fact that 



gsvs- 



gsvs 



< ■T\\gs\\H--'\\vs\\Hi(n) ^ as (5 ^ 0, by ( |3.2[ ), 



we obtain, by (1.3), 



L2(Q2\f^l) 



as (5 ^ 0. 



(3.4) 



Applying Lemma [2| we have 
2 



AVvs • r] 



int 



AVvs • f] 



ext 



< C\\AVvs\\L^[n^\a^){\\AVvs\\L-2{n2\ni) + II div(^Vu5)||i2(f^2\ni) 
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Since div{AVvs) — £g ^eok^T^vs = ^gs in 1^2 \ it follows from (3.2), (3.4), and the 
boundedness of A and S that 



AS/vs • rj 



+ 



This implies 



since esAVvs 



AVvs ■ rj 
■ rj = esAVvs 

ll'^<5ll//i(f^\f^2) - ^ 



cxt 



AS/vs ■ rj 



as (5 — )• 0, 



as (5 ^ 0. 



(3.5) 



• r] on ■ By Lemma jsj we have 
AVvs • rj 



— > as 5 — ^ 0. 



As a consequence, 

11^511^1/2(9^2) ^ as 5 ^ 0. 
Combining (3.4) and (3.7), and using Lemma [T| we have 



P'5||L2(Q2\ni) 



— ;> as (5 — ;> 0. 



We deduce from (3.4) and (3.8) that 

\\vs\\m{n2\n^) ^ as 5 ^ 0. 

This implies 

11^511^1/2(9^1) ^ as (5 ^ 0. 
By Lemma[T| it follows from (3.2) and (3.10) that 



l^5llL2(ni) 



0. 



(3.6) 
(3.7) 
(3.8) 

(3.9) 
(3.10) 

(3.11) 



Multiplying equation ( |3.3| ) by vs and integrating the obtained expression on we have 

-/ 

On the other hand, 



AVvs ■VV5+ / £5 {AVvs, Vu^) - / 



k'^eo'^lvsl'^ 



I 



gsvs- (3.12) 



AVvs ■ rj vs ^ as 6 0, by (3.5) and (3.10); 



^vsr ^ as (5^ 0, by (3.11); 



and 



95VS 



< \\95\\mni)\\v5\\mni) ^ as 5 ^ 0, by ([3^ 



(3.13) 
(3.14) 

(3.15) 
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A combination of (|3.12|), (|3.13|), (|3.14|), and ( |3.15D yields 



From (3.11) and (3.16), we obtain 



as 5 ^ 0. 



(3.16) 



(3.17) 



A combination of (3.6), (3.9), and (3.17) yields 

lb<5|lHi(f^) -J" as (5 : 



we have a contradiction. The proof of (3.1) completes. □ 

3.2 Proof of Theorem H 

We follow the strategy in the proof of Theorem [TJ 
Case 1: / is compatible with the system. The proof is derived from the following steps: 
Step 1: Let v G Hq(Q) be a solution to the equation 

div{eoA\/v) + /c^eoSv = eof. 

We prove that v = NL{k, /). 

Step 2: Define uq := NI{k, /). We prove that uq £ (i^) is a solution to the equation 

div(eo^Vtio) + /c^eo^^^o = £o/- 

Step 3: Prove that {us)o^s<5o is bounded in H^{i^) if 5o is small enough. 
Step 4: Prove that (us) converges weakly in H^{Q,) and strongly in L^(J7) to uq as S 
goes to 0. 

We note immediately here that as a consequence of Steps 1 and 2, NI{k, f) is the 
unique solution in Hq{CI) to the equation (1.14). We will only establish the proof of Steps 
1 and 3. The proof of Steps 2 and 4 is similar to the one of Steps 2 and 4 in Case 1 of the 
proof of Theorem [T] The details are left to the reader. 

Step 1: Assume v £ Hq{Q) is a solution to the equation 



Set 
and 



div(eo^V'u) + fe^eoSf = eof in Q. 



vi = V o F in \ 



1 if G \ 0,2, 
-1 if a; G O4 \ O3. 



(3.18) 
(3.19) 
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It follows from Proposition [l] that vi G H^{^3, \ ^2) H Hj^^^Q^ \ $12) and vi satisfies 

div(eoF*^Vt;i) + k'^ioF^^vi = ioFJ in \ ^2, (3.20) 

and 



= f on and F^AVvi • 7] = AVv • r] 

In the last identity, we use the fact that F*^Vvi • rj = —AVv ■ r] 



on dil.2- 



int 



and AVv ■ rj 



ext 



on 8^2 by Proposition jlj 



AVv ■ r] on dQ2 by the transmission condition on 51^2 • Define 



int 



U = fi — u in \ r22. 
Since F^A = A and = S in \ $^25 it follows that 

div(AVU) + A;2sU = F*/-/ '111^3X^2, 
U = on 5^2, 

AVU-?? = onai72- 



(3.21) 



(3.22) 



Prom the unique continuation principle, see e.g., |22| and and the fact that A G 
0^(0,3 \ $12); we have 

U = [/ in \ ^^2. 
Define the function V2 on Q as follows 

vi o G^^ if X G 0,3, 

V if X G \ Jls- 



V2ix) 



(3.23) 



(3.24) 



Using (1.1) and applying Proposition [T| we have 

div(iV?;2) + k^^V2 = / in \ 8^3 

and, on d^ls, 



(3.25) 



AVv2 ■ rj 



ext 



AVv2 • rj 



AVv ■ T] + F^AVvi ■ T] 



AVv ■ r] 
AVv ■ 7] 



ext 



F^AVvi ■ T] 



ext 



int 



(by (3.20)) 



- F^AV{v + V)-j] (by (3.21)). 

ext int 



It follows from (3.23) that 

AVv2 ■ rj — AVv2 ■ Tj 



ext 



int 



-F^AVU ■ T] 



int 



Since G{x) = x on 8Q, we also obtain, on 8^3 

V2 



- V2 


= V 


- Vl 


= V 


- Vl 


= V 


-{U + v) 


= -u 




ext 


int 


ext 


ext 


ext 


int 


ext 


int 


int 



(3.26) 



(3.27) 
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A combination of (|3.25|), (|3.26|), (|3.27|), and (|1.21l) yields 



^2 = in 17. 

This implies that u satisfies 

' div(AV^;) + k'^T.v = f 
v = W 
AVv ■ T] 



(3.28) 



ext 



int 



AVW -r] 



ext 



in \ ^2, 
on Sils, 
on dQ^,- 



Since in the Sd case ^ G C^(r23 \ 172)7 applying the unique continuation principle, see e.g. 
[22] and [1], we have 

V = y in \ O2. (3.29) 



From (3.24), (3.21), (3.18), it follows that 



V2 in Q. \ O3, 

\] + vi in 173 \ 172, 

vio F in 172 \ 17i, 

V2 o G o F in 17i. 



We deduce from (3.23), (3.28) and (3.29) that 



W in 17 \ 173, 

V in 173 \ 172, 

{V + U)oF inl72\17i, 

WoGoF inl7i, 

Here in the third equality, we also used ( 3.21[ ) and (3.23). Therefore, 

V = NL{k,f) in 17. 



Step 3: Set 



V5 



Us — uq in 17. 
We have, in 17, 

div{esAVvs) + k'^eoT.vs = div{esAVus) — div{esAVuo) + k'^eoT.us — k'^eoT.uo- 
Since dw{esAVuo) + k'^eo'^UQ = eof and div(eo^Vuo) + k'^eoTiUo = eof, it follows that 
diY{esVvs) + k^eoTiVs = div [(eq — £s)AVuo\ in 17. 
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Applying Proposition [2| we obtain 

l|Vw5||L2(f2) < c(||Vuo||l2(q) + ||div(^Vno)||i2(t^2\j^^) + || div(^V'Uo)||L2(t^^) 

<C(||Vno||i2 + ||/|U2(f,)), 
which yields, since us = vs + uq, 

W'^ushnn) < C(||Vno||L2 + \\f\\mn))- 
Since us £ Hq{Q), by the Poincare inequality, it follows that 

\\us\\m{n) < C'(l|Vno||L2 + llfh^n))- 
The proof of Step 3 completes. 

Case 2: / is not compatible with the system. The proof of this case is the same as the one 
of Case 2 in the proof of Theorem [T] The details are left to the reader. □ 

Remark 15. Our analysis (mainly in Steps 1 and 2) would make media whose some 
region contains negative index materials behave like a folded space as noted in IT^ . 

4 Applications 

This section is devoted to some applications of Theorems [l] and [2] for super lenses, 
cloaking, and illusion optics. The use of negative index for super lenses, cloaking, and 
illusion optics are not new as mentioned Section The goal of this section is to present 
a more general setting for these applications whose the rigorous analysis is firmly estab- 
lished in the previous sections, hence to illuminate the general mechanism behind these 
phenomena. 



4.1 Applications for super lenses 



The first immediate application of Theorems [T] and [2] is to create perfect lenses. More 
precisely, consider the settings in Theorems [T] and [2} Assume / is compatible and / = 
in S^s. Since uq = uo G o F on dVti^ one can focus light from d^i to d^^; more precisely 
one has a "perfect lens". 

Here are concrete examples. Let < ri < r2 < rs, m, / > 1, VLi = Br,, 0,2 = B 



= with r3 



Oa 



ro X 



/Ixp and G{x) 



ri 1 ''2 — J-^r2 ) 

x\K 



We recall 



here that Bj. denotes the ball centered at the origin of radius r. We have 



GoF{x) 
Choose m,l > 1 such that 



m{l-l) 



|^|m/-m-Z^ 



ml = m + I. 
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We obtain 

GoF{x) 

It follows that if / is compatible then 



'2 ^ 



(4.1) 



uo{x) = uo(^x) on dBri or mo(x) = ^of ^x) on 5S 



We have a perfect lens. We also have, from (4.1) and (1.7), 

G^F^A{x) -- 

and 



A(rrx/r^) in ^,.3 



(4.2) 



^(rT'x/r^) in Br,. 



Consider the case d = 3, k = 0, A = I in Q\ Br2, 

A = ^rm Br, \ B^^i^^_, , and A = p-^A in Br, \ B 



We have 



A = I = F^A = G^F^A in Br, \ B 



T2- 



Assume / = for /r'l^ ^ < |x| < (this implies V = W). Then / is compatible and 



Uq{X) = Uo 



In particular. 



Uo[X) = UQ 



^x] in B, 



Pendry in [21], and Anantha Ramakrishna and Pendry in |24j (see also [23]) considered 
the case m = I = 2 and / is a dipole. We show here how to obtain similar phenomena for 
a class of transformations F and G. 

Consider the case d = 2, k = 0, and ^ = / in i7 \ il2i A = F~^A in Br, \ Br,, and A 
is constant in f^i. From (4.2), we have 

G^F^A = ^ in Br,. 

Thus the cylindrical Br, would to an outside observer appear magnified by a factor = 
r^/r™~^. This phenomena was first observed by Nicorovici-McPhedran-Milton [T7| in the 
case j4 = / in \ r^i, A = cl in Qi for some constant c > 0, and m = I = 2. 
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4.2 Applications for cloaking and illusion optics 

In this section, we discuss the possibihty to obtain illusion optics and cloaking using 
negative index materials. Given A, T, and A' , S' respectively in \ and il^, we choose 
A in ^2 such that A = F'^A and S = F^^S in O2 \ and A = F~'^G-^A' and 
S = F^^G^^Tj' in Qi. By Theorems [l] and J2j for observers outside Qs, would be 
characterized by the medium A',T,' since A,T, = A',T,' in ^3. Roughly speaking, A,T, 
would appear like A' , S': one obtains illusion optics. By choosing A' , S' = 1, 1 in ^3: one 
cloaks the region \ $12- Here we assume that sources are compatible. The mechanism 
of illusion optics and cloaking described here was respectively proposed in [11] and [l^ in 
the case r^i, fl2, and Q3 F, G, are radial symmetry. These mechanisms have supported 
by simulations and experiments. As far as we know, our paper provides the first rigorous 
analysis for the mechanism of cloaking and illusion optics given in these two papers. 
Nevertheless, we note that in [10], the authors only assumed that A = F~^A and S = 
F~^Ti in Q2\^i; we do not know how to explain this situation using our analysis. The first 
observation on the cloaking using negative index materials (for the 2d-quasi-static case) 
went back to [T7] . Another interesting phenomena related to the mechanism of cloaking, 
we mentioned here, is the cloaking effects associated with anomalous localized resonance 
discussed in a very interesting paper of Milton and Nicorovici [13] (see also [15] , [18] [5] , 
[3j for related works). They showed that in the 2d quasi static regime, with f^i = i?^^, 
Q2 = Br2, and A = I, for a dipole source inside i? 3/2 , 2 \ -6^2, ||i*(5||Hi(n) go^s to infinity; 

r2 /rj V ; 

however us converges outside B 3/2 , 2- This result has been recently extended in [2] where 
they characterized / such that this phenomena happens for arbitrary Qi and 0.2- In this 
paper, we provide a general class of / (/ is compatible with the system) such that this 
phenomena does not take place for a general setting. It would be nice to characterize / 
so that one obtains the cloaking effects associated with anomalous localized resonance i.e. 
1 1 ^5 1 1 /^^{Q) go6s to 00 while as us converges outside a large region. We emphasize here that 
there is a difference in our characterization and the characterization given in [2|. In this 
paper, we characterize sources for which the total energy of the fields remains bounded 
as the loss parameter goes to zero. In {2}, the authors characterize sources for which the 
dissipation energy blows up and the fields remains bounded outside a bounded region as 
the loss parameter goes to zero. 

A Appendix: Scattered case 

In this appendix, we present the corresponding results for scattered waves. For the 
convenience of the reader, we will recall the notation presented in Section [l.l[ Let Q2, 
Q^, and be smooth subsets of M'^ (d = 2, 3) such that they are connected, 

ill CC 0.2 CC O3 cc 0.4. 

We assume that there exist xi G f^i and X2 € and two diffeomorphisms 

F:02\{a;i}^r?4\f^2, G : Q^K^s ^ ^3\{x2} 
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such that there exist diffeomorphism extension of F and G on a neighborhood of and 
d^ls respectively, 

F{n2\^i) = n3\n2, 

F(x) = X on 8^2 and G{x) = x on dQ^, 
and 

G o F : ill — )• is a diffeomorphism if one sets G o F{xi) = X2- 
The geometry is illuminated in Figure 1. 
For 5 >0, es is given as follows. 

'' -l + i5 if \ 1^1 



es{x) :-- 



1 



otherwise , 



We will assume that A and E satisfy (1.3), (1.4), 

A(x) = I and = 1 if |x| is large enough, 

and in the 3d case, we assume in addition that A £ C^{Cl^ \ 0.2)- 

In this appendix, we investigate the behavior of ug G ifj^^^(M"'), as 6 goes to 0, the 
solution to the equation 

dW{esAVus) + k'^eo^us = eof in M"', (Al) 

which satisfies the outgoing condition: 

drUs — ikus = o{r^^^'^^^'^) as r = — )• 00. 

Definition Al. We say that f G L^(M^) (d = 2 or 3) with compact support is compatible 
with the system if there exist U G H^{i^3 \ Q2) o-nd V G H'^{Q.z \ O2) such that 

' dw{AVU) + k'^m = FJ - f in n3\^2, 

U = on 8^2, 

AVU-i] = on 8^2, 

and 

' div{Avv) + k^j:v = f in n3\n2, 

V = W on d^s, 

ext 

AW-r] = AVW-r] on 8^3. 

cxt 

Here W G r\jiyoH^{Bji \ 8^1.3) is the unique outgoing solution to the system 

' div(iVW^) + k^tW = f in R*^ \ (9^3 , 
[W] = -U on 8^3 , 

[AVW • 7]] = -AVU • rj on 8^3 , 
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If f is compatible with the system, then NIs{k,f) is defined as follows 



V if x€ 9.3X^2, 

{u + v)oF if xen2\^i, 
WoGoF ifxeni. 



(A2) 



NIs{kJ) = 
We are ready to state 

Theorem Al. Let 6 > 0, d = 2, or 3, / G L'^{W^) with compact support, and us G 
H^^^{M.'^) be the unique outgoing solution to the equation (Al). Assume that F^:A = A and 
F^T, = Tj in ils \ r22, and in the 3d case, assume in addition that A G C^{Q3 \ ^^2)- We 
have 

a) Case 1: f is compatible with the system. The sequence {us) converges weakly in 
Hl^{W^) and strongly in if^^(M^) to uq G HI^{W^), the unique outgoing solution to 
the equation 

div(eoVtio) + k^SolluQ = e^f in M''. 



as 5 — )• 0. Moreover, uq = NLs{k,f) where NLs{k,f) is given in (A2). 
b) Case 2: f is not compatible with the system. We have 

lim 11^511^1^1(5^) = +00 if 0,3, supp/ CC Br. 



Assume in addition that f = in Q3 and G*-F*74 = yl in \ Then U = and 
V = W m^l3\Q2, and W G //^^^^(M*^) is the outgoing solution to the equation 



dw{AVW) + k'^J:W = / in 



We have the following 



Corollary A3. Assume G^F^A = F^A = ^ in \ O2, ^ G C^{Q.3 \ 9.2), and f = in 



^3. Then f is compatible and (us), the solutions to the equations (Al), converges weakly 
in Hj-^^(K'^) and strongly in L^^^(M'^) to uo the unique outgoing solution in Hj^^^CR.'^) to the 
equation 

div(eoAV'Uo) + k'^eoT,uo = eof in W^, 

as (5 — )• 0. Moreover 

uo = W in M'^\173, 
where W G Hj-^^{M.'^) is the unique outgoing solution to the equation 

div(iVH^) + k'^tW = f in M^. 

A cousin of Theorem |3] in which the loss parameter involves in the two terms is as 
follows: 
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Theorem A2. Let 5 > d = 2, or 3, f ^ L^(M'^) with compact support, and us G 
Hj-^^{]R.'^) be the unique outgoing solution to the equation 

div{esAVus) + k^esT^us = f in R'^. 

Assume that F^:A = A and = S in ^3 \ il2; (ind in the 3d case, assume in addition 
that A G C'^{n3\n2)- We have 

a) Case 1: f is compatible with the system. The sequence {us) converges weakly in 
Hl^iR"^) and strongly in L^^iR'^) to uq G HI^{R'^), the unique outgoing solution to 
the equation 

div(eoVtio) + k'^eo'SuQ = eof in R'^. 



as 5 ^ 0. Moreover, uq = NLs{k,f) where NLs{k,f) is given in (A2). 
b) Case 2: f is not compatible with the system. We have 

lim 11^511^1^1(5^) = +00 if 0,3, supp/ CC Br. 
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